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Abstract. We study the zero temperature Casimir energy and fermion number 
for Dirac fields in a 2+1-dimensional Minkowski space-time, in the presence of a 
uniform magnetic field perpendicular to the spatial manifold. Then, we go to the 
finite-temperature problem, with a chemical potential, introduced as a uniform 
zero component of the gauge potential. By performing a Lorentz boost, we obtain 
Hall's conductivity in the case of crossed electric and magnetic fields. 
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1. Introduction 

The quantization of the Hall conductivity ^ is a remarkable quantum phenomenon, 
which occurs in two-dimensional electron systems, at low temperatures and strong 
perpendicular magnetic fields. Most proposed explanations for this phenomenon [2] 
rely on Schroedinger's one-particle theory and the introduction of a filling fraction of 
the Landau levels, which must be assumed to be integer to reproduce the observed 
behavior. 

It is the aim of this paper to show that, in the context of relativistic field theory, 
such behavior naturally arises, as a consequence of the spin-statistics theorem. 

In section |21 we present the theory of Dirac fields in 2 + 1 Minkowski space-time, 
interacting with a magnetic background field perpendicular to the spatial plane, and 
evaluate the vaccum expectation values of the energy and fermion density. Section 
121 contains the generalities of the same theory in Euclidean 3-dimcnsional space, and 
presents the eigenvalues of the corresponding Dirac operator. From such eigenvalues, 
the partition function is evaluated in section 01 Section contains the resulting 
free energy and mean particle density at finite temperature. Finally, in section 
we perform an adequate Lorentz boost in order to consider the problem of fermions 
interaction with crossed electric and magnetic fields, and obtain the Hall conductivity. 
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2. Zero-temperature problem 

We study a 2 + 1-dimensional theory of Dirac fields, in the presence of a uniform 
background magnetic field perpendicular to the spatial plane. We choose the metric 
(—,+,+), natural units K = c = 1, and adopt the following representation for the 
Dirac matrices: 7^ = ias, 7^ = 02 ^-^Im ~ 

The Hamiltonian can be determined from the solutions of the Dirac equation 
{if)— eA)^> = 0, where — e is the negative charge of the electron. In the Landau gauge 
A = (0,0, Bx) , with B > 0. Thus, after setting ^{t,x, y) = e~ iEt tf>(x,y), we get the 
Hamiltonian H = i<7id x — io2^y + (Xi&Bx. 

In order to solve the eigenvalue problem for this Hamiltonian, we take 

This leads to the following system of first order equations 
(id x - ik - ieBx)xk = Eip k 

(id x + ik + ieBx)tp k — E\k ■ (2) 
After imposing that the eigenfunctions be well-behaved for x — > ±00, we find two 
types of solutions to our problem ||: 

1. A definite-chirality zero mode (Eq = 0) 

2. A set of eigenfunctions corresponding to the symmetric spectrum 
E n = ±^/2neB , n = 1, 00. 

In all cases, the eigenfunctions can be written in terms of Hermite polynomials, 
and all eigenvalues exhibit the well known Landau degeneracy per unit area: 

A,=g ; (3) 

The vacuum expectation value of the energy per unit area, defined through a 
zeta function regularization (see, for example, 3 and references therein), is given by 

Ec = -^Y,E n ^\En\- S \ ■ 

J s=— 1 

In the present case, we have (a is an arbitrary parameter with mass dimension, 
introduced to render the complex powers dimensionless) 



a 



1 



= -A L V2eBQ R (-- ) . (4) 



s=-l 

Always in the zeta-function regularization framework, the fermion number is 0j 

n(b)=~^(yi \ e ^ s - e i^n +^0, 

\£„>0 £„<0 / J s=0 

where Nq is the contribution coming from zero modes. 

In our case, the nonvanishing spectrum is symmetric. So, only the zero mode, 
which is charge self-conjugate, contributes. This gives as a result 

N(B) = ±^. (5) 

Or, equivalently, for the vacuum expectation value of the charge density 

f(B)= T e^. (6) 



|| For the other non-equivalent representation of the gamma matrices in 2+1 dimensions, the spectrum 
doesn't change. 
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3. The theory at finite temperature with chemical potential 

In order to study the effect of temperature, we go to Euclidean space, with the metric 
(+, +, +) . To this end, we take the Euclidean gamma matrices to be 70 = i-ff.j = —03, 
71 = 7^ = (72, 72 = — 0i- We will follow |S| in introducing the chemical potential 
as an imaginary Aq = —i^ in Euclidean space. Thus, the partition function in the 
grand-canonical ensemble is given by 

\nZ = In det(ip- eA) . (7) 

In order to evaluate it in the zeta regularization approach [H], we first determine 
the eigenfunctions, and the corresponding eigenvalues, of the Dirac operator, in the 
same gauge used in the previous section, i.e, we solve 

[-ia 3 (d T + n) + ia 2 d x + oi{id v - eBi)]f = , (8) 

To satisfy antiperiodic boundary conditions in the t direction, we propose 

e iX t T e iky n 

* k>l (T,x,y) = Mx), with A, = (21 + 1)-, (9) 

y/2np p 

where (3 — ^ is the inverse temperature. 

After doing so, and writing ibh Ax) — ( ^ }k ' l ^f : \ ] we have for each fc, L 

V XkA x ) J 

(d x - k - eBx)xk,i =(u-\i)(fk,i 

(-d x - k - eBx)ipk,i = {u> + Xi)xk,i , (10) 

where A/ = A; — ifi. 

There are two types of eigenvalues 
1. uji = A/, with I — —00, 00. 



2. u>i n = ±y Xf + 2neB, with n = 1, 00, I = —00, 00. 

In all cases, the degeneracy per unit area is again given by Al in equation @. 



4. Evaluation of the partition function at finite temperature and chemical 
potential 



The partition function, in the zeta regularization scheme [3], is given by 

d 



log Z 



ds 



(11) 



s=0 



As in the previous section, a is a parameter with mass dimension, introduced to render 
the ^-function dimensionless. 

We must consider two contributions to log Z, respectively coming from eigenvalues 
of type 1 and 2 in the previous section, i.e., 

Ai(/x) = ~ Cx (*,/*) 

s=0 



ds 

d_ 

ds 



l=-c 



ap a 



(12) 
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and 



A 2 (/x,B)= --r 
as 



d 
ds 



( 2 (s,n,B) 



s=0 



(l+(-l)-)Ai 



OO 

E 

71=1 

' = — OO 



2neB 



+ (21 + 1) 



a/3 



(13) 



In the rest of this section, we sketch the main steps in the analytic extension 
of both zeta functions and in the calculation of their s-derivatives (for a detailed 
presentation, see [7]). The contribution Ai(fi) can be evaluated at once for the whole 
/i-range. The analytic extension of Ci( s ) A*) can be achieved as follows (for a similar 
calculation, see [S]) 



Ci(«,M) = A L 



2tt 

a/3 



1 ifi/3 
*' 2 ~ ~2n~ 



E 

1=0 



(14) 



Now, in order to write the second term as a Hurwitz zeta, we must relate the 
eigenvalues with negative real part to those with positive one without, in so doing, 
going through zeros in the argument of the power. Otherwise stated, we must select 
a cut in the complex u> plane UJ. This requirement determines a definite value of 
(— l)~ s , i.e., (— l)~ s = e l7rs *s™(p) s . Taking this into account, we finally have 



Ci(s,M) = A L ( 

1 pa 



1 



CH ( flj ±_^) +e *^>. Cff ( 8) ± + ^ 



2tt 



1 



ifi/3 



2tt 



(15) 



From this last expression, the contribution Ai(/x) to logZ can be obtained. It is 
given by 



AiM = - Ai 



1 



i/i/3 



2tt 



A L { log I 2 cosh ( ^ 



i/i/3 



2tt 



2tt 



(16) 



The analytic extension of £2(3, /*t, -B) requires a separate consideration of different 
H ranges. We study in detail two of these ranges 

4.1. fi 2 < 2eB 



( 2 ( S ,^B) = (l + (-l)- s )A L 



OO 



2neB 



n = 1 L 

/ = — 00 



(a + i)^-<£ 

ap a 



(17) 



Making use of the Mellin transform, this can be written as 



( 2 {s,^,B) 



(1 + (-!)-«) A, 

r(f) 

-2t { IT 



E 



_-, Jo 



n=l 
ill 



dtt^~ l e 



4nt 
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where we have used the definition of the Jacobi theta function 

To proceed, we use the inversion formula for the Jacobi function, 
03(2, x) = -^e^^^03 (#:, i), and perform the integration over t, thus getting 



C2 (s,fi,B) = 



(l + (-l)- s )A L af3s 

4V^r(^) 



s - 1 



2eS 



s- 1 



4 B-D 



nJ=l 



l 2 q 4 /? 2 
8« eJ B 



cosh 



(/x/3i) A'^_i [y2neBf3 2 t- 



(19) 



From this expression, the contribution A2 to the partition function can be readily 
obtained, since the factor accompanying s is finite at s = 0. After using that 
K_i[x) = \p^e~ x ' 1 and performing the resulting sum over I, we obtain 

r 



A 2 (fx,B) = A L (3 

+ ^log(l + e" 2Vra/3 + 2cosh( A i/3) e - v ^ 3/3 ) 



(20) 



Finally, adding the contributions given by equations <|16f) and (|20|l we get, for the 
partition function in the range /1 2 < 2eB 



\ogZ = A L < log 2 cosh 



M/3 
2 



OO 

+ £ lo § (l + e- 2v ^ /3 + 2 cosh {nP)e-^*P) 



(21) 



2eS < n 2 < 4eB 
As before, we have 

00 

C 2 ( S ,A 1 ,S) = (l + (-l)- s )A L J2 

n = 1 L 

l = —00 



2neB 



(2/ + !)^-^ 

ap a 



(22) 



However, in this range of /Li, the contribution to the zeta function due to n = 1 
must be analytically extended in a different way. In fact, the expression cannot be 
written in terms of a unique Mellin transform, since its real part is not always positive 
(note, in connection with this that, for n = 1, eq. (|19|) diverges). Instead, it can be 
written as a product of two Mellin transforms 



,n=l ( r>\ (l + Mr S ) 

C 2 [s, Ml B) - 



^ poo 

A L J2 d ^f-i e -[( 2 '+i)f-v+"V2Z 



B I 



(23) 



Relativistic Landau problem 



or, after changing variables according to t' = t — z; z' = t + z, performing one of the 
integrals, and the sum over I 



dzz 2 Js-i (V2eBz 



-fi . 



(24) 



'sinh(f) 

Now, the integral in this expression diverges at z = 0. In order to isolate such 
divergence, we add and subtract the first term in the series expansion of the Bessel 
function, thus getting the following two pieces 

(l + (-l)- s )^F s 



4a- s r(§ + 1) 



dz , 



J^i{V2e~Bz) - 



sinh(^) 



and 

fn=U q „ m_ (l + (-l)- fl )Vi 



dzz 8 - 1 



(25) 



(26) 



sinh(f ) 

The contribution of equation 1(25(1 to the partition function can be easily evaluated 
by noticing that the factor multiplying s is finite at s — 0. This gives 

-2M/3 



A^ 1 1 ) ( A1 ,B) = A L {log(l 



2cosh(V2^B/3)e 



>l/3 



|//| i 2l«H;|2cosh(^) 



(27) 



In order to get the contribution coming from 1)26(1 . the integral can be evaluated 
for 5fts > 1, which gives 

(1 + (-l)- s )r( s )0F(«/3) s A L 



(2 7 r) s 2 s - 1 r(§)r(^±i) 
Ch(s, x(i )) + Ch(s, tt(i + — )) 

Z 7T Z TT 



(28) 



where £y(s,2;) is the Hurwitz zeta function. Its contribution to the partition function 
can now be evaluated by using that ( H (0, ±(1 - lJ ^-)) + Ch(0, |(1 + = and the 
well known value of — s =o(h (s, x) ^H], to obtain 

'%)). (29) 

Summing up the contributions in equations (|16|l . H27|) and (|29|l . as well as the 
contribution coming from n > 2, evaluated as in the previous subsection, one gets for 
the partition function 



A^ 1 , (fi, B) = 2A L log(2 cosh (^)) . 



K»-Z = A z . Jlo,, ( zco.l, 



log(l 



,-2H/9 



2cosh(\/2^B/3)e- |Al|/3 ) + /V2^B (Cfl(-^J -1 
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+ J2 l °S (l + e- 2 ^^? + 2 cosh ( M /3) e -^^) . (30) 

At first sight, this result looks different from the one corresponding to /x 2 < 2eB 
(equation J2J). However, it is easy to see that both expressions coincide. The 
advantage of using expression (|30|l for this range of fi is that the zero-temperature 
limit is explicitly isolated from finite-temperature corrections. 



5. Free energy and particle number 



From equations l|21(l and (|30|l , the free energy per unit area (F = — 4 log Z) can be 
obtained \ It is given by 



(31) 

, oo. In the 



F(n,B,0)= -A £ |^ log (2 cosh I -if + V2eB( R \-~ 

+ i log (l + e - 2 v^S/3 + 2 cosh ( M(5 ) e -^B^ , 

71— 1 J 

Moreover, the free energy is continuous at [x 2 — 2neB, n = 0, . 
low-temperature limit one has 

F(^ 2 < 2e£) ^^oo -A L V2e~B<; R (- 1 

which coincides with the Casimir energy obtained in sectional even for /x ^ 0, but in 
this range, i.e., for \i less than the first Landau level, if positive, or greater than minus 
the first Landau level, if negative. On the other hand, 

F(2eB < i? < AeB) -A L jV^B Ur (-^j - lj + • 

The mean particle density can be obtained as N = -g^-logZ. For nonzero 
temperature and arbitrary fi (not coinciding with an energy level" 1 ") one has 



N(n,B,0)=A L t- 



1 r tanh(^) -tignM 



2sinh (n/3)e- V2l ™ s P 



(32) 



^ l + e -2V2^B(3 + 2 cos h ( Ai/ 3) e -v / 2^/3 } 

It is interesting to note that, for fi = 0, one has iV(/i = 0) = ±4^, which shows 
that, in the absence of chemical potential, the fermion number obtained in equation 
|5| remains unaltered as the temperature grows. 

On the other hand, for nonvanishing fi, the low-temperature limit differs, 
depending on the grange considered 

N(2eBn < ii 2 < 2eB(n + 1)) — > /3 ^ 00 nA L sign(fi) , 

r o 

where n 



2eB 



% Consistently with the footnote in section [3] all the results in this section are independent from the 
representation of the gamma matrices chosen. 

+ Note that, for instance, if n = 2neB, the series in equation 1191 converges only conditionally, and 
its term- by-term derivative leads to a divergent series. 
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This result is nothing but the expected one for particles with the statistic of 
fermions, since relativistic field theory naturally leads to the spin-statistics theorem. 
At zero temperature, fi is nothing but the Fermi energy; for example, for /i > 0, as fi 
grows past a Landau level, such level becomes entirely filled. 



6. Final comments 



From the previous result, the mean value of the particle density at zero temperature 
can be obtained. After recovering units, one has 

f(2ec 2 hBn < y 2 < 2eBc 2 h(n + 1)) = - sign(fi) , 

the other two components of the current density tri-vector being equal to zero in the 

absence of an electric field. 

Now, the zero-temperature limit of the same tri-vector in the presence of crossed 

homogeneous electric (F 1 ) and magnetic (B') fields can retrieved, for F' < cB' , by 

performing a Lorentz boost with absolute value of the velocity ^7. Suppose, for 

definiteness, that the homogeneous electric field points along the positive y axis. 

Then, the velocity of the Lorentz boost must point along the negative x-axis, and 

the transformation gives as a result 

.,0 ~nce 2 B' . , /X -ne 2 F' , .,„ 
J = 7 sign{n) , j = sign(^) , j = . 

As a consequence, the quantized zero-temperature Hall conductivity is 
— ne 
~~h 

Finally, we mention that the more realistic case of massive fermions is at present 
under study [TT], 



v xy = — — sign{p) . 



Acknowledgments 

We thank Paola Giacconi and Roberto Soldati for useful discussions. 

This work was partially supported by Universidad Nacional de La Plata, under Grant 

11/X381. 



References 

[1] K. v. Klitzing. The quantized Hall effect. Rev. Mod. Phys., 58:519, 1986. 

[2] T. Chakraborty and P. Pietilainen. Quantum Hall effects: Integral and fractional. Springer- 
Verlag, Berlin, 1995. 

[3] E. Elizalde. Ten physical applications of spectral zeta functions. Springer- Verlag, Berlin, 1995. 
[4] A.J. Niemi and G.W. Semenoff. Fermion Number Fractionization in Quantum Field Theory. 

Phys. Reports, 135:99, 1986. 
[5] A. Actor. Chemical Potentials in Gauge Theories. Phys. Lett., B157:53, 1985. 
[6] J.S. Dowker and R. Critchley. Effective Lagrangian and energy momentum tensor in de Sitter 

space. Phys. Rev., D13:3224, 1976. 
[7] C.G. Beneventano and E.M. Santangelo. in preparation. 2005. 

[8] C.G. Beneventano and E.M. Santangelo. Finite temperature properties of the Dirac operator 
under local boundary conditions. Jour. Phys. A: Math. Gen., 37:9261, 2004. 

[9] G. Cognola, E. Elizalde and S. Zerbini. Dirac Functional Determinants in Terms of the Eta 
Invariant and the Noncommutative Residue. Commun. Math. Phys. 237:507, 2003. 
[10] L.S. Gradshteyn and L.M. Ryzhik. Table of Integrals, Series and Products. Academic Press, 
USA, 2000. 

[11] C.G. Beneventano, P. Giacconi, E.M. Santangelo and R. Soldati. Work in progress. 2005. 



